Practice of teacher's example¶
In [10]:
library(tidyverse)
─ Attaching packages ──────────────────── tidyverse 1.2.1 ─ ✔ ggplot2 2.2.1 ✔ purrr 0.2.4 ✔ tibble 1.4.2 ✔ dplyr 0.7.4 ✔ tidyr 0.8.0 ✔ stringr 1.3.0 ✔ readr 1.1.1 ✔ forcats 0.3.0 ─ Conflicts ───────────────────── tidyverse_conflicts() ─ ✖ dplyr::filter() masks stats::filter() ✖ dplyr::lag() masks stats::lag()
In [11]:
diamonds
| carat | cut | color | clarity | depth | table | price | x | y | z |
|---|---|---|---|---|---|---|---|---|---|
| 0.23 | Ideal | E | SI2 | 61.5 | 55 | 326 | 3.95 | 3.98 | 2.43 |
| 0.21 | Premium | E | SI1 | 59.8 | 61 | 326 | 3.89 | 3.84 | 2.31 |
| 0.23 | Good | E | VS1 | 56.9 | 65 | 327 | 4.05 | 4.07 | 2.31 |
| 0.29 | Premium | I | VS2 | 62.4 | 58 | 334 | 4.20 | 4.23 | 2.63 |
| 0.31 | Good | J | SI2 | 63.3 | 58 | 335 | 4.34 | 4.35 | 2.75 |
| 0.24 | Very Good | J | VVS2 | 62.8 | 57 | 336 | 3.94 | 3.96 | 2.48 |
| 0.24 | Very Good | I | VVS1 | 62.3 | 57 | 336 | 3.95 | 3.98 | 2.47 |
| 0.26 | Very Good | H | SI1 | 61.9 | 55 | 337 | 4.07 | 4.11 | 2.53 |
| 0.22 | Fair | E | VS2 | 65.1 | 61 | 337 | 3.87 | 3.78 | 2.49 |
| 0.23 | Very Good | H | VS1 | 59.4 | 61 | 338 | 4.00 | 4.05 | 2.39 |
| 0.30 | Good | J | SI1 | 64.0 | 55 | 339 | 4.25 | 4.28 | 2.73 |
| 0.23 | Ideal | J | VS1 | 62.8 | 56 | 340 | 3.93 | 3.90 | 2.46 |
| 0.22 | Premium | F | SI1 | 60.4 | 61 | 342 | 3.88 | 3.84 | 2.33 |
| 0.31 | Ideal | J | SI2 | 62.2 | 54 | 344 | 4.35 | 4.37 | 2.71 |
| 0.20 | Premium | E | SI2 | 60.2 | 62 | 345 | 3.79 | 3.75 | 2.27 |
| 0.32 | Premium | E | I1 | 60.9 | 58 | 345 | 4.38 | 4.42 | 2.68 |
| 0.30 | Ideal | I | SI2 | 62.0 | 54 | 348 | 4.31 | 4.34 | 2.68 |
| 0.30 | Good | J | SI1 | 63.4 | 54 | 351 | 4.23 | 4.29 | 2.70 |
| 0.30 | Good | J | SI1 | 63.8 | 56 | 351 | 4.23 | 4.26 | 2.71 |
| 0.30 | Very Good | J | SI1 | 62.7 | 59 | 351 | 4.21 | 4.27 | 2.66 |
| 0.30 | Good | I | SI2 | 63.3 | 56 | 351 | 4.26 | 4.30 | 2.71 |
| 0.23 | Very Good | E | VS2 | 63.8 | 55 | 352 | 3.85 | 3.92 | 2.48 |
| 0.23 | Very Good | H | VS1 | 61.0 | 57 | 353 | 3.94 | 3.96 | 2.41 |
| 0.31 | Very Good | J | SI1 | 59.4 | 62 | 353 | 4.39 | 4.43 | 2.62 |
| 0.31 | Very Good | J | SI1 | 58.1 | 62 | 353 | 4.44 | 4.47 | 2.59 |
| 0.23 | Very Good | G | VVS2 | 60.4 | 58 | 354 | 3.97 | 4.01 | 2.41 |
| 0.24 | Premium | I | VS1 | 62.5 | 57 | 355 | 3.97 | 3.94 | 2.47 |
| 0.30 | Very Good | J | VS2 | 62.2 | 57 | 357 | 4.28 | 4.30 | 2.67 |
| 0.23 | Very Good | D | VS2 | 60.5 | 61 | 357 | 3.96 | 3.97 | 2.40 |
| 0.23 | Very Good | F | VS1 | 60.9 | 57 | 357 | 3.96 | 3.99 | 2.42 |
| ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
| 0.70 | Premium | E | SI1 | 60.5 | 58 | 2753 | 5.74 | 5.77 | 3.48 |
| 0.57 | Premium | E | IF | 59.8 | 60 | 2753 | 5.43 | 5.38 | 3.23 |
| 0.61 | Premium | F | VVS1 | 61.8 | 59 | 2753 | 5.48 | 5.40 | 3.36 |
| 0.80 | Good | G | VS2 | 64.2 | 58 | 2753 | 5.84 | 5.81 | 3.74 |
| 0.84 | Good | I | VS1 | 63.7 | 59 | 2753 | 5.94 | 5.90 | 3.77 |
| 0.77 | Ideal | E | SI2 | 62.1 | 56 | 2753 | 5.84 | 5.86 | 3.63 |
| 0.74 | Good | D | SI1 | 63.1 | 59 | 2753 | 5.71 | 5.74 | 3.61 |
| 0.90 | Very Good | J | SI1 | 63.2 | 60 | 2753 | 6.12 | 6.09 | 3.86 |
| 0.76 | Premium | I | VS1 | 59.3 | 62 | 2753 | 5.93 | 5.85 | 3.49 |
| 0.76 | Ideal | I | VVS1 | 62.2 | 55 | 2753 | 5.89 | 5.87 | 3.66 |
| 0.70 | Very Good | E | VS2 | 62.4 | 60 | 2755 | 5.57 | 5.61 | 3.49 |
| 0.70 | Very Good | E | VS2 | 62.8 | 60 | 2755 | 5.59 | 5.65 | 3.53 |
| 0.70 | Very Good | D | VS1 | 63.1 | 59 | 2755 | 5.67 | 5.58 | 3.55 |
| 0.73 | Ideal | I | VS2 | 61.3 | 56 | 2756 | 5.80 | 5.84 | 3.57 |
| 0.73 | Ideal | I | VS2 | 61.6 | 55 | 2756 | 5.82 | 5.84 | 3.59 |
| 0.79 | Ideal | I | SI1 | 61.6 | 56 | 2756 | 5.95 | 5.97 | 3.67 |
| 0.71 | Ideal | E | SI1 | 61.9 | 56 | 2756 | 5.71 | 5.73 | 3.54 |
| 0.79 | Good | F | SI1 | 58.1 | 59 | 2756 | 6.06 | 6.13 | 3.54 |
| 0.79 | Premium | E | SI2 | 61.4 | 58 | 2756 | 6.03 | 5.96 | 3.68 |
| 0.71 | Ideal | G | VS1 | 61.4 | 56 | 2756 | 5.76 | 5.73 | 3.53 |
| 0.71 | Premium | E | SI1 | 60.5 | 55 | 2756 | 5.79 | 5.74 | 3.49 |
| 0.71 | Premium | F | SI1 | 59.8 | 62 | 2756 | 5.74 | 5.73 | 3.43 |
| 0.70 | Very Good | E | VS2 | 60.5 | 59 | 2757 | 5.71 | 5.76 | 3.47 |
| 0.70 | Very Good | E | VS2 | 61.2 | 59 | 2757 | 5.69 | 5.72 | 3.49 |
| 0.72 | Premium | D | SI1 | 62.7 | 59 | 2757 | 5.69 | 5.73 | 3.58 |
| 0.72 | Ideal | D | SI1 | 60.8 | 57 | 2757 | 5.75 | 5.76 | 3.50 |
| 0.72 | Good | D | SI1 | 63.1 | 55 | 2757 | 5.69 | 5.75 | 3.61 |
| 0.70 | Very Good | D | SI1 | 62.8 | 60 | 2757 | 5.66 | 5.68 | 3.56 |
| 0.86 | Premium | H | SI2 | 61.0 | 58 | 2757 | 6.15 | 6.12 | 3.74 |
| 0.75 | Ideal | D | SI2 | 62.2 | 55 | 2757 | 5.83 | 5.87 | 3.64 |
Single Variable : Category¶
In [12]:
ggplot(data = diamonds, aes(x = cut)) + geom_bar(fill = "yellow", color = "black")
Single Variable : Continuous¶
In [13]:
ggplot( data = diamonds, aes(x = price)) + geom_histogram(binwidth = 10)
Two Variable : Continuous vs Continuous¶
In [14]:
ggplot( data = diamonds, aes(x = table, y = depth)) + geom_point()
Two variable and Two variable vs. Continuous¶
In [47]:
ggplot(data = diamonds, aes(x = clarity, y = depth)) + geom_boxplot()
Many variables¶
In [15]:
library(ggplot2)
library(GGally)
library(scales)
library(memisc)
Attaching package: ‘GGally’
The following object is masked from ‘package:dplyr’:
nasa
Attaching package: ‘scales’
The following object is masked from ‘package:purrr’:
discard
The following object is masked from ‘package:readr’:
col_factor
Loading required package: lattice
Loading required package: MASS
Attaching package: ‘MASS’
The following object is masked from ‘package:dplyr’:
select
Attaching package: ‘memisc’
The following object is masked from ‘package:scales’:
percent
The following objects are masked from ‘package:dplyr’:
collect, recode, rename
The following objects are masked from ‘package:stats’:
contr.sum, contr.treatment, contrasts
The following object is masked from ‘package:base’:
as.array
ggpairs 用法
In [19]:
# 從數據集獲取10000個樣本數據進行分析
set.seed(20022012)
# lists that may contain the variables 'continuous', 'combo', 'discrete', and 'na'
diamond_samp <- diamonds[sample(1:length(diamonds$price), 10000), ]
ggpairs(diamond_samp,lower= list(continuous = wrap("points", shape = I('.'))),
upper = list(combo = wrap("box", outlier.shape = I('.'))))
`stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`. `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
geom_abline_with_ggploty¶
In [20]:
library(plotly)
packageVersion('plotly')
Attaching package: ‘plotly’
The following objects are masked from ‘package:memisc’:
rename, style
The following object is masked from ‘package:MASS’:
select
The following object is masked from ‘package:ggplot2’:
last_plot
The following object is masked from ‘package:stats’:
filter
The following object is masked from ‘package:graphics’:
layout
[1] ‘4.7.1’
In [21]:
set.seed(1234)
rnorm( n, mean = k, sd = g)
n : 隨機變數個數; mean : 平均數; sd : 標準差
In [24]:
dat <- data.frame(cond = factor(rep(c("A", "B"), each=200)), rating = c(rnorm(200)))
In [25]:
dat
| cond | rating |
|---|---|
| A | -1.22681509 |
| A | 0.03615287 |
| A | -0.42139311 |
| A | -0.89936441 |
| A | 0.41744132 |
| A | 0.15344474 |
| A | 1.46328305 |
| A | -1.12150250 |
| A | -0.51778808 |
| A | -0.07494709 |
| A | -1.40779008 |
| A | -0.28470587 |
| A | -0.70817118 |
| A | -2.14763900 |
| A | -0.28383716 |
| A | -0.53407216 |
| A | 1.13301019 |
| A | -0.60406892 |
| A | 0.55751160 |
| A | 0.14262929 |
| A | -1.23686021 |
| A | 0.37414397 |
| A | -0.10608612 |
| A | 0.16606199 |
| A | 0.65803485 |
| A | 0.11396292 |
| A | 1.58515916 |
| A | 0.11555344 |
| A | -1.25399904 |
| A | 0.74186478 |
| ⋮ | ⋮ |
| B | -0.22276341 |
| B | -1.77146803 |
| B | 0.28742225 |
| B | -0.64119176 |
| B | -0.92493764 |
| B | -0.77997747 |
| B | 0.44450337 |
| B | -2.16246900 |
| B | -0.42551090 |
| B | -1.01371766 |
| B | -0.96211903 |
| B | 0.47463463 |
| B | 0.35175427 |
| B | -0.10702036 |
| B | 0.53751726 |
| B | -0.60210490 |
| B | -2.02329765 |
| B | 1.15442128 |
| B | -0.23069300 |
| B | -0.36508877 |
| B | 0.15751793 |
| B | 0.76593222 |
| B | 0.42608697 |
| B | 0.92726540 |
| B | 0.04613459 |
| B | -0.50428564 |
| B | -0.30044941 |
| B | 0.82291598 |
| B | -0.84926880 |
| B | -0.95361010 |
In [26]:
p <- ggplot(dat, aes(x = rating)) +
geom_histogram(binwidth = .25, fill = "grey", color = "black") +
geom_vline(aes(xintercept = mean(rating, na.rm=T)), color = "red", linetype = "dashed", size = .5)
## na.rm=T ==> remove na value
In [27]:
p <- ggplotly(p)
We recommend that you use the dev version of ggplot2 with `ggplotly()`
Install it with: `devtools::install_github('hadley/ggplot2')`
In [28]:
p
require is used inside functions, as it outputs a warning and continues if the package is not found
library will throw an error if not found
In [3]:
require(graphics)
In [2]:
typeof(UCBAdmissions)
'double'
In [7]:
UCBAdmissions[]
<NA>
In [21]:
UCBAdmissions
, , Dept = A
Gender
Admit Male Female
Admitted 512 89
Rejected 313 19
, , Dept = B
Gender
Admit Male Female
Admitted 353 17
Rejected 207 8
, , Dept = C
Gender
Admit Male Female
Admitted 120 202
Rejected 205 391
, , Dept = D
Gender
Admit Male Female
Admitted 138 131
Rejected 279 244
, , Dept = E
Gender
Admit Male Female
Admitted 53 94
Rejected 138 299
, , Dept = F
Gender
Admit Male Female
Admitted 22 24
Rejected 351 317
In [22]:
apply(UCBAdmissions, c(1, 2), sum)
| Male | Female | |
|---|---|---|
| Admitted | 1198 | 557 |
| Rejected | 1493 | 1278 |
apply 函數家族¶
apply 是Ｒ中核心數據處理的函數。可以實現數據的循環、分組、過濾、類型操作。¶
- apply() : 參數list, data.frame, array 返回vector, matrix
- tapply() : 參數vector 返回vector
- vapply() : 參數list, data.frame 返回vector, matrix
- lapply() : 參數list, data.frame 返回list
- sapply() : 參數list, data.frame 返回vector, matrix
- mapply() : 參數vector, 返回vector, matrix
- rapply() : 參數list, 返回list
- eapply() : 參數environment, 返回list
In [8]:
## Data aggregated over departments
# iterate through data
apply(UCBAdmissions, c(1, 2), sum)
mosaicplot(apply(UCBAdmissions, c(1, 2), sum),
main = "Student admissions at UC Berkeley")
## Data for individual departments
opar <- par(mfrow = c(2, 3), oma = c(0, 0, 2, 0))
for(i in 1:6)
mosaicplot(UCBAdmissions[,,i],
xlab = "Admit", ylab = "Sex",
main = paste("Department", LETTERS[i]))
mtext(expression(bold("Student admissions at UC Berkeley")),
outer = TRUE, cex = 1.5)
par(opar)
| Male | Female | |
|---|---|---|
| Admitted | 1198 | 557 |
| Rejected | 1493 | 1278 |
In [24]:
require(stats)
require(graphics) # for mosaicplot
In [25]:
summary(esoph)
agegp alcgp tobgp ncases ncontrols 25-34:15 0-39g/day:23 0-9g/day:24 Min. : 0.000 Min. : 1.00 35-44:15 40-79 :23 10-19 :24 1st Qu.: 0.000 1st Qu.: 3.00 45-54:16 80-119 :21 20-29 :20 Median : 1.000 Median : 6.00 55-64:16 120+ :21 30+ :20 Mean : 2.273 Mean :11.08 65-74:15 3rd Qu.: 4.000 3rd Qu.:14.00 75+ :11 Max. :17.000 Max. :60.00
In [26]:
model1 <- glm(cbind(ncases, ncontrols) ~ agegp + tobgp * alcgp,
data = esoph, family = binomial())
In [27]:
model1
Call: glm(formula = cbind(ncases, ncontrols) ~ agegp + tobgp * alcgp,
family = binomial(), data = esoph)
Coefficients:
(Intercept) agegp.L agegp.Q agegp.C
-1.75985 2.99646 -1.35008 0.13436
agegp^4 agegp^5 tobgp.L tobgp.Q
0.07098 -0.21347 0.63846 0.02922
tobgp.C alcgp.L alcgp.Q alcgp.C
0.15607 1.37077 -0.14913 0.22823
tobgp.L:alcgp.L tobgp.Q:alcgp.L tobgp.C:alcgp.L tobgp.L:alcgp.Q
-0.70426 0.12225 -0.29187 0.12948
tobgp.Q:alcgp.Q tobgp.C:alcgp.Q tobgp.L:alcgp.C tobgp.Q:alcgp.C
-0.44527 -0.05205 -0.16118 0.04843
tobgp.C:alcgp.C
-0.13905
Degrees of Freedom: 87 Total (i.e. Null); 67 Residual
Null Deviance: 227.2
Residual Deviance: 47.48 AIC: 237
In [28]:
anova(model1)
| Df | Deviance | Resid. Df | Resid. Dev | |
|---|---|---|---|---|
| NULL | NA | NA | 87 | 227.24056 |
| agegp | 5 | 88.12808 | 82 | 139.11248 |
| tobgp | 3 | 19.08475 | 79 | 120.02773 |
| alcgp | 3 | 66.05445 | 76 | 53.97328 |
| tobgp:alcgp | 9 | 6.48949 | 67 | 47.48379 |
In [29]:
model2 <- glm(cbind(ncases, ncontrols) ~ agegp + unclass(tobgp)
+ unclass(alcgp),
data = esoph, family = binomial())
In [30]:
model2
Call: glm(formula = cbind(ncases, ncontrols) ~ agegp + unclass(tobgp) +
unclass(alcgp), family = binomial(), data = esoph)
Coefficients:
(Intercept) agegp.L agegp.Q agegp.C agegp^4
-4.01097 2.96113 -1.33735 0.15292 0.06668
agegp^5 unclass(tobgp) unclass(alcgp)
-0.20288 0.26162 0.65308
Degrees of Freedom: 87 Total (i.e. Null); 80 Residual
Null Deviance: 227.2
Residual Deviance: 59.28 AIC: 222.8
In [31]:
summary(model2)
Call:
glm(formula = cbind(ncases, ncontrols) ~ agegp + unclass(tobgp) +
unclass(alcgp), family = binomial(), data = esoph)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.7628 -0.6426 -0.2709 0.3043 2.0421
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -4.01097 0.31224 -12.846 < 2e-16 ***
agegp.L 2.96113 0.65092 4.549 5.39e-06 ***
agegp.Q -1.33735 0.58918 -2.270 0.02322 *
agegp.C 0.15292 0.44792 0.341 0.73281
agegp^4 0.06668 0.30776 0.217 0.82848
agegp^5 -0.20288 0.19523 -1.039 0.29872
unclass(tobgp) 0.26162 0.08198 3.191 0.00142 **
unclass(alcgp) 0.65308 0.08452 7.727 1.10e-14 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 227.241 on 87 degrees of freedom
Residual deviance: 59.277 on 80 degrees of freedom
AIC: 222.76
Number of Fisher Scoring iterations: 6
In [32]:
ttt <- table(esoph$agegp, esoph$alcgp, esoph$tobgp)
o <- with(esoph, order(tobgp, alcgp, agegp))
ttt[ttt == 1] <- esoph$ncases[o]
tt1 <- table(esoph$agegp, esoph$alcgp, esoph$tobgp)
tt1[tt1 == 1] <- esoph$ncontrols[o]
tt <- array(c(ttt, tt1), c(dim(ttt),2),
c(dimnames(ttt), list(c("Cancer", "control"))))
mosaicplot(tt, main = "esoph data set", color = TRUE)
In [37]:
library(ggplot2)
library(maptools)
library(rgeos)
library(Cairo)
library(ggmap)
library(scales)
library(RColorBrewer)
library(mapproj)
Loading required package: maps
In [42]:
map <- get_map(location = 'Taiwan', zoom = 7)
ggmap(map)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=Taiwan&zoom=7&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Taiwan&sensor=false
In [43]:
# use language to set the word on the map
map <- get_map(location = 'Taiwan', zoom = 7,
language = "zh-TW")
ggmap(map)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=Taiwan&zoom=7&size=640x640&scale=2&maptype=terrain&language=zh-TW&sensor=false Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Taiwan&sensor=false
In [64]:
# to set accurate location
map <- get_map(location = c(lon = 121.533937, lat = 25.03933), zoom = 12, language = "zh-TW")
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=25.03933,121.533937&zoom=12&size=640x640&scale=2&maptype=terrain&language=zh-TW&sensor=false
In [65]:
ggmap(map)
In [59]:
# Different type
map <- get_map(location = c(lon = 121.533937, lat = 25.03933), zoom = 13, language = "zh-TW", maptype = "roadmap")
ggmap(map)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=25.03933,121.533937&zoom=13&size=640x640&scale=2&maptype=roadmap&language=zh-TW&sensor=false
In [66]:
# Great to show data on the map
map <- get_map(location = c(lon = 121.533937, lat = 25.03933), zoom = 13, language = "zh-TW", maptype = "toner-lite")
ggmap(map)
maptype = "toner-lite" is only available with source = "stamen". resetting to source = "stamen"... Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=25.03933,121.533937&zoom=13&size=640x640&scale=2&maptype=terrain&sensor=false Map from URL : http://tile.stamen.com/toner-lite/13/6860/3505.png Map from URL : http://tile.stamen.com/toner-lite/13/6861/3505.png Map from URL : http://tile.stamen.com/toner-lite/13/6862/3505.png Map from URL : http://tile.stamen.com/toner-lite/13/6860/3506.png Map from URL : http://tile.stamen.com/toner-lite/13/6861/3506.png Map from URL : http://tile.stamen.com/toner-lite/13/6862/3506.png Map from URL : http://tile.stamen.com/toner-lite/13/6860/3507.png Map from URL : http://tile.stamen.com/toner-lite/13/6861/3507.png Map from URL : http://tile.stamen.com/toner-lite/13/6862/3507.png Map from URL : http://tile.stamen.com/toner-lite/13/6860/3508.png Map from URL : http://tile.stamen.com/toner-lite/13/6861/3508.png Map from URL : http://tile.stamen.com/toner-lite/13/6862/3508.png
In [1]:
data <- read.csv("UV_20180319014346.csv")
In [2]:
data
| SiteName | UVI | PublishAgency | County | WGS84Lon | WGS84Lat | PublishTime |
|---|---|---|---|---|---|---|
| 花蓮 | 0.00 | 中央氣象局 | 花蓮縣 | 121,36,48 | 23,58,30 | 2018-03-19 01:00 |
| 馬祖 | 0.00 | 中央氣象局 | 連江縣 | 119,55,24 | 26,10,09 | 2018-03-19 01:00 |
| 高雄 | 0.00 | 中央氣象局 | 高雄市 | 120,18,57 | 22,33,58 | 2018-03-19 01:00 |
| 玉山 | 0.00 | 中央氣象局 | 南投縣 | 120,57,34 | 23,29,15 | 2018-03-19 01:00 |
| 臺南 | 0.00 | 中央氣象局 | 臺南市 | 120,12,17 | 22,59,36 | 2018-03-19 01:00 |
| 新竹 | 0.00 | 中央氣象局 | 新竹縣 | 121,00,51 | 24,49,40 | 2018-03-19 01:00 |
| 鞍部 | 0.00 | 中央氣象局 | 臺北市 | 121,31,47 | 25,10,57 | 2018-03-19 01:00 |
| 恆春 | 0.00 | 中央氣象局 | 屏東縣 | 120,44,47 | 22,00,14 | 2018-03-19 01:00 |
| 臺北 | 0.00 | 中央氣象局 | 臺北市 | 121,30,53 | 25,02,16 | 2018-03-19 01:00 |
| 成功 | 0.00 | 中央氣象局 | 臺東縣 | 121,22,24 | 23,05,51 | 2018-03-19 01:00 |
| 基隆 | 0.00 | 中央氣象局 | 基隆市 | 121,44,26 | 25,07,60 | 2018-03-19 01:00 |
| 新屋 | 0.00 | 中央氣象局 | 桃園市 | 121,02,51 | 25,00,24 | 2018-03-19 01:00 |
| 蘭嶼 | 0.00 | 中央氣象局 | 臺東縣 | 121,33,30 | 22,02,13 | 2018-03-19 01:00 |
| 臺東 | 0.00 | 中央氣象局 | 臺東縣 | 121,09,17 | 22,45,08 | 2018-03-19 01:00 |
| 日月潭 | 0.00 | 中央氣象局 | 南投縣 | 120,54,29 | 23,52,53 | 2018-03-19 01:00 |
| 金門 | 0.00 | 中央氣象局 | 金門縣 | 118,17,21 | 24,24,26 | 2018-03-19 01:00 |
| 宜蘭 | 0.00 | 中央氣象局 | 宜蘭縣 | 121,45,24 | 24,45,50 | 2018-03-19 01:00 |
| 澎湖 | 0.06 | 中央氣象局 | 澎湖縣 | 119,33,47 | 23,33,56 | 2018-03-19 01:00 |
| 臺中 | 0.00 | 中央氣象局 | 臺中市 | 120,41,03 | 24,08,45 | 2018-03-19 01:00 |
| 嘉義 | 0.00 | 中央氣象局 | 嘉義市 | 120,25,58 | 23,29,45 | 2018-03-19 01:00 |
| 塔塔加 | NA | 環境保護署 | 嘉義縣 | 120,52,50.06 | 23,28,14.19 | 2018-03-19 01:00 |
| 阿里山 | NA | 環境保護署 | 嘉義縣 | 120,48,05.02 | 23,30,30.82 | 2018-03-18 23:00 |
| 屏東 | NA | 環境保護署 | 屏東縣 | 120,29,16.92 | 22,40,23.09 | 2018-03-19 00:00 |
| 橋頭 | NA | 環境保護署 | 高雄市 | 120,18,20.48 | 22,45,27.02 | 2018-03-19 00:00 |
| 新營 | NA | 環境保護署 | 臺南市 | 120,19,2.10 | 23,18,20.28 | 2018-03-19 00:00 |
| 朴子 | NA | 環境保護署 | 嘉義縣 | 120,14,50.46 | 23,27,55.11 | 2018-03-19 00:00 |
| 斗六 | NA | 環境保護署 | 雲林縣 | 120,32,41.98 | 23,42,42.67 | 2018-03-19 00:00 |
| 南投 | NA | 環境保護署 | 南投縣 | 120,41,7.1 | 23,54,46.8 | 2018-03-19 00:00 |
| 彰化 | NA | 環境保護署 | 彰化縣 | 120,32,29.47 | 24,03,57.60 | 2018-03-19 00:00 |
| 沙鹿 | NA | 環境保護署 | 臺中市 | 120,34,7.66 | 24,13,32.26 | 2018-03-19 00:00 |
| 苗栗 | NA | 環境保護署 | 苗栗縣 | 120,49,12.72 | 24,33,54.97 | 2018-03-19 00:00 |
| 桃園 | NA | 環境保護署 | 桃園市 | 121,19,11.87 | 24,59,41.24 | 2018-03-19 00:00 |
| 板橋 | NA | 環境保護署 | 新北市 | 121,27,31.2 | 25,00,46.7 | 2018-03-19 00:00 |
| 淡水 | NA | 環境保護署 | 新北市 | 121,26,57.26 | 25,09,52.20 | 2018-03-19 00:00 |
In [7]:
lon.deg <- sapply((strsplit(as.character(data$WGS84Lon), ",")), as.numeric)
data$lon <- lon.deg[1, ] + lon.deg[2, ]/60 + lon.deg[3, ]/3600
lat.deg <- sapply((strsplit(as.character(data$WGS84Lat), ",")), as.numeric)
data$lat <- lat.deg[1, ] + lat.deg[2, ]/60 + lat.deg[3, ]/3600
In [8]:
lon.deg
| 121 | 119 | 120 | 120 | 120 | 121 | 121 | 120 | 121 | 121 | ⋯ | 120.0 | 120.00 | 120.00 | 120.0 | 120.00 | 120.00 | 120.00 | 121.00 | 121.0 | 121.00 |
| 36 | 55 | 18 | 57 | 12 | 0 | 31 | 44 | 30 | 22 | ⋯ | 19.0 | 14.00 | 32.00 | 41.0 | 32.00 | 34.00 | 49.00 | 19.00 | 27.0 | 26.00 |
| 48 | 24 | 57 | 34 | 17 | 51 | 47 | 47 | 53 | 24 | ⋯ | 2.1 | 50.46 | 41.98 | 7.1 | 29.47 | 7.66 | 12.72 | 11.87 | 31.2 | 57.26 |
In [6]:
data$lat
- 23.975
- 26.1691666666667
- 22.5661111111111
- 23.4875
- 22.9933333333333
- 24.8277777777778
- 25.1825
- 22.0038888888889
- 25.0377777777778
- 23.0975
- 25.1333333333333
- 25.0066666666667
- 22.0369444444444
- 22.7522222222222
- 23.8813888888889
- 24.4072222222222
- 24.7638888888889
- 23.5655555555556
- 24.1458333333333
- 23.4958333333333
- 23.4706083333333
- 23.5085611111111
- 22.6730805555556
- 22.7575055555556
- 23.3056333333333
- 23.4653083333333
- 23.7118527777778
- 23.913
- 24.066
- 24.2256277777778
- 24.5652694444444
- 24.9947888888889
- 25.0129722222222
- 25.1645
In [4]:
map <- get_map(location = 'Taiwan', zoom = 7)
Error in get_map(location = "Taiwan", zoom = 7): 沒有這個函數 "get_map" Traceback:
In [81]:
ggmap(map) + geom_point(aes(x = lon, y = lat, size = UVI), data = data)
Warning message: “Removed 14 rows containing missing values (geom_point).”
In [83]:
ggmap(map) +
geom_point(aes(x = lon, y = lat, size = UVI), data = data) +
facet_grid( ~ PublishAgency)
Warning message: “Removed 14 rows containing missing values (geom_point).”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <92>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <a2>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <83>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bf>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <9d>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b7>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <e7>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <bd>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '環境保護署' in 'mbcsToSbcs': dot substituted for <b2>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e4>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b8>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ad>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a4>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <ae>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e6>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b0>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a3>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e8>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <a1>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <e5>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <b1>”Warning message in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, : “conversion failure on '中央氣象局' in 'mbcsToSbcs': dot substituted for <80>”
In [84]:
d <- function(x=-95.36, y=29.76, n,r,a){
round(data.frame(
lon = jitter(rep(x,n), amount = a),
lat = jitter(rep(y,n), amount = a)
), digits = r)
}
df <- d(n = 50,r = 3,a = .3)
map <- get_googlemap(markers = df, path = df,, scale = 2)
ggmap(map)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=29.763284,-95.363271&zoom=10&size=640x640&scale=2&maptype=terrain&markers=29.625,-95.119%7C29.461,-95.463%7C29.895,-95.55%7C29.764,-95.241%7C29.851,-95.519%7C29.828,-95.414%7C29.629,-95.224%7C29.976,-95.392%7C29.624,-95.222%7C29.862,-95.429%7C29.612,-95.483%7C29.881,-95.502%7C29.754,-95.443%7C29.586,-95.539%7C29.47,-95.54%7C29.789,-95.441%7C29.518,-95.112%7C29.614,-95.493%7C30.033,-95.558%7C29.592,-95.581%7C29.719,-95.456%7C29.588,-95.512%7C29.653,-95.293%7C29.976,-95.117%7C29.485,-95.459%7C29.848,-95.623%7C29.808,-95.497%7C29.761,-95.197%7C29.736,-95.292%7C29.756,-95.15%7C29.516,-95.623%7C29.851,-95.177%7C29.925,-95.163%7C29.828,-95.208%7C29.593,-95.322%7C29.853,-95.363%7C29.633,-95.543%7C29.472,-95.453%7C29.659,-95.495%7C29.974,-95.276%7C30.041,-95.572%7C29.801,-95.148%7C29.878,-95.277%7C29.759,-95.183%7C29.619,-95.333%7C29.86,-95.643%7C29.695,-95.155%7C29.563,-95.43%7C29.466,-95.096%7C29.874,-95.44&path=29.625,-95.119%7C29.461,-95.463%7C29.895,-95.55%7C29.764,-95.241%7C29.851,-95.519%7C29.828,-95.414%7C29.629,-95.224%7C29.976,-95.392%7C29.624,-95.222%7C29.862,-95.429%7C29.612,-95.483%7C29.881,-95.502%7C29.754,-95.443%7C29.586,-95.539%7C29.47,-95.54%7C29.789,-95.441%7C29.518,-95.112%7C29.614,-95.493%7C30.033,-95.558%7C29.592,-95.581%7C29.719,-95.456%7C29.588,-95.512%7C29.653,-95.293%7C29.976,-95.117%7C29.485,-95.459%7C29.848,-95.623%7C29.808,-95.497%7C29.761,-95.197%7C29.736,-95.292%7C29.756,-95.15%7C29.516,-95.623%7C29.851,-95.177%7C29.925,-95.163%7C29.828,-95.208%7C29.593,-95.322%7C29.853,-95.363%7C29.633,-95.543%7C29.472,-95.453%7C29.659,-95.495%7C29.974,-95.276%7C30.041,-95.572%7C29.801,-95.148%7C29.878,-95.277%7C29.759,-95.183%7C29.619,-95.333%7C29.86,-95.643%7C29.695,-95.155%7C29.563,-95.43%7C29.466,-95.096%7C29.874,-95.44&sensor=false
In [85]:
mu <- c(-95.3632715, 29.7632836)
nDataSets <- sample(4:10,1)
chkpts <- NULL
for(k in 1:nDataSets){
a <- rnorm(2); b <- rnorm(2);
si <- 1/3000 * (outer(a,a) + outer(b,b))
chkpts <- rbind(chkpts,
cbind(MASS::mvrnorm(rpois(1,50), jitter(mu, .01), si), k))
}
chkpts <- data.frame(chkpts)
names(chkpts) <- c("lon", "lat","class")
chkpts$class <- factor(chkpts$class)
qplot(lon, lat, data = chkpts, colour = class)
In [86]:
ggmap(get_map(maptype = "satellite"), extent = "device") +
stat_density2d(aes(x = lon, y = lat, colour = class), data = chkpts, bins = 5)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=29.763284,-95.363271&zoom=10&size=640x640&scale=2&maptype=satellite&language=en-EN&sensor=false Warning message: “`panel.margin` is deprecated. Please use `panel.spacing` property instead”
In [89]:
# only violent crimes
violent_crimes <- subset(crime,
offense != "auto theft" &
offense != "theft" &
offense != "burglary"
)
# rank violent crimes
violent_crimes$offense <-
factor(violent_crimes$offense,
levels = c("robbery", "aggravated assault",
"rape", "murder")
)
# restrict to downtown
violent_crimes <- subset(violent_crimes,
-95.39681 <= lon & lon <= -95.34188 &
29.73631 <= lat & lat <= 29.78400
)
In [90]:
library(grid)
theme_set(theme_bw(16))
HoustonMap <- qmap("houston", zoom = 14, color = "bw")
# a contour plot
HoustonMap +
stat_density2d(aes(x = lon, y = lat, colour = offense),
size = 3, bins = 2, alpha = 3/4, data = violent_crimes) +
scale_colour_discrete("Offense", labels = c("Robery","Aggravated Assault","Rape","Murder")) +
theme(
legend.text = element_text(size = 15, vjust = .5),
legend.title = element_text(size = 15,face="bold"),
legend.key.size = unit(1.8,"lines")
)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=houston&zoom=14&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=houston&sensor=false Warning message: “`panel.margin` is deprecated. Please use `panel.spacing` property instead”Warning message: “Removed 11 rows containing non-finite values (stat_density2d).”
In [91]:
# 二維的 histogram
HoustonMap +
stat_bin2d(aes(x = lon, y = lat, colour = offense, fill = offense),
size = .5, bins = 30, alpha = 2/4, data = violent_crimes) +
scale_colour_discrete("Offense",
labels = c("Robery","Aggravated Assault","Rape","Murder"),
guide = FALSE) +
scale_fill_discrete("Offense", labels = c("Robery","Aggravated Assault","Rape","Murder")) +
theme(
legend.text = element_text(size = 15, vjust = .5),
legend.title = element_text(size = 15,face="bold"),
legend.key.size = unit(1.8,"lines")
)
Warning message: “Removed 11 rows containing non-finite values (stat_bin2d).”Warning message: “Removed 8 rows containing missing values (geom_tile).”
In [96]:
HoustonMap +
stat_density2d(aes(x = lon, y = lat, fill = ..level.., alpha = ..level..),
size = 2, bins = 4, data = violent_crimes, geom = "polygon") +
scale_fill_gradient("Violent\nCrime\nDensity") +
scale_alpha(range = c(.4, .75), guide = FALSE) +
guides(fill = guide_colorbar(barwidth = 1.5, barheight = 10))
Warning message: “Removed 11 rows containing non-finite values (stat_density2d).”
In [98]:
houston <- get_map("houston", zoom = 14)
overlay <- stat_density2d(aes(x = lon, y = lat, fill = ..level.., alpha = ..level..), bins = 4, geom = "polygon", data = violent_crimes)
HoustonMap +
stat_density2d(aes(x = lon, y = lat, fill = ..level.., alpha = ..level..),
bins = 4, geom = "polygon", data = violent_crimes) +
scale_fill_gradient("Violent\nCrime\nDensity") +
scale_alpha(range = c(.4, .75), guide = FALSE) +
guides(fill = guide_colorbar(barwidth = 1.5, barheight = 10)) +
inset(
grob = ggplotGrob(ggplot() + overlay +
scale_fill_gradient("Violent\nCrime\nDensity") +
scale_alpha(range = c(.4, .75), guide = FALSE) +
theme_inset()
),
xmin = attr(houston,"bb")$ll.lon +
(7/10) * (attr(houston,"bb")$ur.lon - attr(houston,"bb")$ll.lon),
xmax = Inf,
ymin = -Inf,
ymax = attr(houston,"bb")$ll.lat +
(3/10) * (attr(houston,"bb")$ur.lat - attr(houston,"bb")$ll.lat)
)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=houston&zoom=14&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=houston&sensor=false
Warning message:
“Removed 11 rows containing non-finite values (stat_density2d).”Warning message in min(x):
“min 中沒有無漏失的引數; 回傳 Inf”Warning message in max(x):
“max 中沒有無漏失的引數；回傳 -Inf”Warning message in min(x):
“min 中沒有無漏失的引數; 回傳 Inf”Warning message in max(x):
“max 中沒有無漏失的引數；回傳 -Inf”ERROR while rich displaying an object: Error in grid.Call.graphics(C_setviewport, vp, TRUE): non-finite location and/or size for viewport
Traceback:
1. FUN(X[[i]], ...)
2. tryCatch(withCallingHandlers({
. rpr <- mime2repr[[mime]](obj)
. if (is.null(rpr))
. return(NULL)
. prepare_content(is.raw(rpr), rpr)
. }, error = error_handler), error = outer_handler)
3. tryCatchList(expr, classes, parentenv, handlers)
4. tryCatchOne(expr, names, parentenv, handlers[[1L]])
5. doTryCatch(return(expr), name, parentenv, handler)
6. withCallingHandlers({
. rpr <- mime2repr[[mime]](obj)
. if (is.null(rpr))
. return(NULL)
. prepare_content(is.raw(rpr), rpr)
. }, error = error_handler)
7. mime2repr[[mime]](obj)
8. repr_text.default(obj)
9. paste(capture.output(print(obj)), collapse = "\n")
10. capture.output(print(obj))
11. evalVis(expr)
12. withVisible(eval(expr, pf))
13. eval(expr, pf)
14. eval(expr, pf)
15. print(obj)
16. print.ggplot(obj)
17. grid.draw(gtable)
18. grid.draw.gTree(gtable)
19. recordGraphics(drawGTree(x), list(x = x), getNamespace("grid"))
20. drawGTree(x)
21. grid.draw(x$children[[i]], recording = FALSE)
22. grid.draw.gTree(x$children[[i]], recording = FALSE)
23. recordGraphics(drawGTree(x), list(x = x), getNamespace("grid"))
24. drawGTree(x)
25. grid.draw(x$children[[i]], recording = FALSE)
26. grid.draw.gTree(x$children[[i]], recording = FALSE)
27. recordGraphics(drawGTree(x), list(x = x), getNamespace("grid"))
28. drawGTree(x)
29. preDraw(x)
30. preDraw.gTree(x)
31. pushvpgp(x)
32. pushgrobvp(x$vp)
33. pushgrobvp.viewport(x$vp)
34. pushViewport(vp, recording = FALSE)
35. lapply(vps, push.vp, recording)
36. FUN(X[[i]], ...)
37. push.vp.vpStack(X[[i]], ...)
38. lapply(vp, push.vp, recording)
39. FUN(X[[i]], ...)
40. push.vp.viewport(X[[i]], ...)
41. grid.Call.graphics(C_setviewport, vp, TRUE)
In [95]:
df <- data.frame(
x = rnorm(10*100, -95.36258, .05),
y = rnorm(10*100, 29.76196, .05),
year = rep(paste("year",format(1:10)), each = 100)
)
for(k in 0:9){
df$x[1:100 + 100*k] <- df$x[1:100 + 100*k] + sqrt(.05)*cos(2*pi*k/10)
df$y[1:100 + 100*k] <- df$y[1:100 + 100*k] + sqrt(.05)*sin(2*pi*k/10)
}
ggmap(get_map(),
base_layer = ggplot(aes(x = x, y = y), data = df)) +
stat_density2d(aes(fill = ..level.., alpha = ..level..),
bins = 4, geom = "polygon") +
scale_fill_gradient2(low = "white", mid = "orange", high = "red", midpoint = 10) +
scale_alpha(range = c(.2, .75), guide = FALSE) +
facet_wrap(~ year)
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=29.763284,-95.363271&zoom=10&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false
